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I) 

Abstract: We study a system of a complex charged scalar coupled to a Reissner-Nordstrom 

^\j , black hole in 3+1 dimensional anti-de Sitter spacetime, neglecting back-reaction. With 

^ \ suitable boundary conditions, the cases of a neutral and purely electric black hole have been 

\^ \ studied in various limits and were shown to yield key elements of superconductivity in the 

^T ' dual 2+1 dimensional field theory, forming a condensate below a critical temperature. By 

adding magnetic charge to the black hole, we immerse the superconductor into an external 
f~^ ■ magnetic field. We show that a family of condensates can form and we examine their structure. 

OO . For finite magnetic field, they are localized in one dimension with a profile that is exactly 

solvable, since it maps to the quantum harmonic oscillator. As the magnetic field increases. 



the condensate shrinks in size, which is reminiscent of the Meissner effect. 
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1. Introduction 

Since the early days of AdS/CFT[Q, |2|, ^ (see ref. [Q] for a review) it has been tempting to 
consider applying holographic duality to the study of important strongly coupled phenomena 
in condensed matter systems, and superconductivty has been high on the list^. Since finite 
temperature in these duals typically implies the presence of a black hole, and since supercon- 
ductivity requires a condensate to form below a certain critical temperature, the existence of a 
holographically dual background would seem to require a circumvention of various statements 
of no-hair theorems (which go back to Wheeler [^] — for a review, see ref.Q). Generically, 
the black holes would need to have some kind of scalar hair in order to be dual to a super- 
conductor (the scalar's asymptotic value would be the condensate vacuum expectation value 
(vev)). 

In a series of studiesjl^, ^, ^], Gubser has presented a case for just the right kind of 
no-scalar-hair theorem evasion to allow for a superconductor's dual to exist. The statement 
seems to be that there do exist solutions that allow for a condensing scalar to be coupled 
to the black hole if the charge on the black hole is large enough. The scalar couples to (at 
least) a U{1) under which the black hole is charged, and its condensation breaks the gauge 
symmetry spontaneously, giving a mass to the gauge field. In particular, if the effective mass 



^The fact that certain black holes and branes are known to exhibit a sort of Meissner effect H, o, M at zero 
temperature has always added to the motivation, although it has not been clear how exactly this could be 
connected to a dual superconductivity. 



in the bulk of the scalar is negative enough, the scalar field develops a non-trivial vev at the 
boundary, giving the gauge field a non-zero mass. 

Studying such solutions is hard to do, since the full equations are coupled and non-linear, 
and so numerical methods, and a number of limits, have been employed in order to extract 
the key physics. Gubser has studied[^ the case of non-Abelian Reissner-Nordstrom black 
holes condensing, and in a simpler model that seems to capture some of the essentials in a 
limit, Hartnoll et.al.,[|l^] have studied a neutral black hole with a charged scalar and Maxwell 
sector that do not back react on the geometry. The latter authors have explored (with the 
aid of that simplifying limit) some of the phenomenology of the condensate as a function of 
temperature and shown that it maps rather well (where the limit can be trusted) to familiar 
features of superconductivity in the dual 2+1 dimensional theory. 

Emboldened by these studies, we explored the case of adding an external magnetic field 
to the system, to see how the condensate behaves^. We have a fully back-reacted electrically 
and magnetically charged Reissner-Nordstrom black hole, and a charged scalar whose back- 
reaction we neglect in our computations. Since the scalar does not back-react, we cannot 
hope to see all of the signature physics of a superconductor in the presence of magnetism, as 
the superconductor is not able to repel the background magnetic field. Instead, we find that 
the condensate generically adjusts itself so as to fill only a strip of finite width in the plane, 
thereby reducing the total magnetic field that threads it. Remarkably, we can solve exactly 
for the profile that it adopts, and we find that, as the magnetic field approaches infinity, the 
condensate shrinks to zero size. 

2. The Background 

We begin by introducing a charged, complex scalar field into the four dimensional Einstein- 
Maxwell action with a negative cosmological constant^: 

S = -l^ [d'^xV^iR+-^ + L'^(-]F^-\d^-igA^i'\^-Vi\^\)]] . (2.1) 



This action contains a term proportional to A^j^A^^^^. This term contributes negatively to 
the effective mass of the charged scalar since the charged black hole will source At- It is 
exactly this term that allows (but does not guarantee) a non-trivial vev for the scalar field 
to form. When a vev is formed, by the usual Higgs-Anderson mechanism, the gauge field 
develops a mass term proportional to Afj^A'^{^^). In the limit where the scalar field ^ does 
not back-react on the geometry, the solution for the background geometry we take is that of 
the dyonic black hole [|l5|] : 



z^ z^ f [z] 



As we were preparing this manuscript, a paper on the same subject (ref. Il4]) appeared on the Ar^iv 
^We are using the mostly positive signature convention. 



F = 2ha^dx f\dy + 2qadz A dt , 
f^,) = l+^h' + q^),^-^l + h' + q^),^ = (l-,)^,2 + , + l-^h^ + q')z^) .(2.2) 

In the coordinate system used in equation (p.2|), z is a dimensionless radial coordinate scaled 
so that the event horizon of the black hole is located at Zh = 1 and the AdS boundary 
is at z ^ 0. The parameters a, h, and q are related to the mass, magnetic charge, and 
electric charge of the black hole respectively, but only a is dimensionful, with dimension of 
inverse length. These quantities are in turn related to the temperature, external magnetic 
field, and charge density of the charged adjoint matter in the dual field theory. The only 
other dimensionful parameter in the solution is L, related to the AdS radius. The Hawking 



temperature is given by the usual Gibbons-Hawking calculus[16|: 



T='- = ^{3-h^-q^). (2.3) 

Note that in order for the temperature to remain positive, (/i^ + q"^) < 3. Saturating this 
inequality corresponds to the extremal, zero-temperature case. In order to determine the 
effect of the magnetic and electric charges of the black hole, we choose a particular form for 
the gauge field A, such that F = dA: 

A = 2ha^xdy + 2qa (z - I) dt . (2.4) 

We have explicitly added the pure gauge term —2qadt in order to have A regular at the event 
horizon jl^ . The Ay term provides a constant magnetic field B = F^y , and this is interpreted 
as corresponding to an external magnetic field in the (2+l)-dimensional system [p^]: 

B = 2ha^ . (2.5) 

The At term has two terms, a constant term and a term that goes to zero at the boundary. 
The constant term is interpreted as the chemical potential (for an analogue of R-charge; see 
e.g., refs[^, ^), and the second term can be related to the conjugate dual charge density of 
the theory via: 

1 "Don-shell ^ 2 

V^Mt(z = 0) ~~K^"" 
where V is the volume of the two-dimensional spatial part of the field theory. 

3. The Scalar Field 

3.1 Review 

Let us review the results of ref. [Of to better clarify the relationship to our present work. In 
that paper, the background is neutral, so both the electric and magnetic charge of the dyonic 
black hole have been set to zero. Instead, the Maxwell-scalar sector is decoupled from the 
gravity sector by sending the coupling g ^ oo. In order to see this, we must first rescale 



J- w'-'on— shell ^2 fn n\ 



A^ — > A^/g and ^ -^ ^ /q- The Maxwell-scalar sector then has an overall 5^^, which when 
sent to infinity, decouples it from the gravity sector. In this analysis, the potential is taken 
to be: 

y(l^l) = -2^*/l2 . (3.1) 

Therefore, one can now study the Maxwell-scalar theory in the black hole background with 
Lagrangian: 

C = --F^ - \d^ - iA^^ + 2^^/L^ (3.2) 

The equation of motion for the fields ^ and A^ are: 

^d^ {^/^G^''' {d^-^ - iA^-^)) + ^^- iG^"'A^ (d^^ - iA^^) = , (3.3) 

^0^ ( V^G^'^G^'^i^Aa) - G^'' {i (^5^^ - 9^^^) + 2A^^^) = , (3.4) 



and that of ^ is simply the complex conjugate of equation ( |3.3D . We take the ansatz: 

^ = ^(z) = <&(z)/L, At = At{z) = aAt{z), (3.5) 

where ^ and At are dimensionless fields. It is then consistent to take the phase of ^ to be 
constant. All other fields are set to zero. Under this ansatz, the equations of motion simplify 
to: 

J ^ \ Q ,T, , -'^ ,T, ;r2 , ^ ,f, n £:i2 T ^ ,t,2 



di^ +[j--) d.^ + p'^Af + ^^ = , dtAt - ^^M, = . (3.6) 

Without presenting the details of the analysis (see ref. |l^), we show in figure |l| the results 
of the variation of an order parameter as the temperature changes. The onset of supercon- 
ductivity occurs for T < T^,. The critical temperature Tc is proportional to the square root 
of the charge density. 

3.2 Pert ur bat ive limit 

We now consider the scalar field as a perturbation about the dyonic black hole background. 
In this analysis, the Maxwell-scalar sector is not isolated from the gravity sector, since the 
Maxwell field has back-reacted on the gravitational background. We use the same potential 



considered in equation (3J), which corresponds to choosing m^L^ = —2 for the scalar field. 
Before proceeding with our analysis, we would like to emphasize the relationship between 
this work and that of ref. [^], which we reviewed in the previous section. We work in the 
limit where the scalar does not backreact on the Maxwell fields, which should correspond 
approximately to taking Af ^ <I> in ref. |l^. From equation (3.5), we see that this limit can 



be accomplished by taking aL -^ 00. In this limit, the charge density diverges, and hence the 
limit corresponds to taking T/T^, -^ 0, i. e. the left most end of the curve in figure 111 This 
argument is further established by the fact that, in this regime. At in the coupled equations 



studied in ref. |13| and reviewed earlier behaves almost linearly. In the dyonic black hole 



background. At is linear. This suggests that the analysis we propose captures the physics at 
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Figure 1: The vev of the A = 2 operator as a function of the temperature for the isolated Maxwell- 
scalar sector studied in ref. |0] . Tc is proportional to the square root of the charge density. Note that 
we use a different normalization, which accounts for the difference of a factor of V2 with ref. hM . 



T/Tc ^^ in ref. |13]. Therefore, for the physics we uncover, we are well below the critical 
temperature. 

The equation of motion for the scalar field is given by: 



1 



-G 



-.du 



GG^"" (d,^ - igA,^)) + ;^^ - igG'^'^A^ (9,* - igA,^) = . (3.7) 



The equation of motion for ^ is simply the complex conjugate of equation (3.7). Using the 
fact that we only have At and Ay, and the only dependence is on the coordinates x and z, 
we consider an ansatz of the form ^ = ^ {x,z). The equation of motion simplifies to: 



^d, (V^G^'d,^) + G^^a^^ + ^^ - Gyyg^Al (x) ^ - G^^g'^Al (z) ^ = . (3. 
—G ^ ^ ^ 



This equation and its complex conjugate are purely real. Therefore, the equations of motion 
imply that the phase of ^ is constant, and so without loss of generality, we take ^ to be real. 
We assume a separable form for ^: 



^! = X{x)Z{z) 
which further simplifies the equation of motion to: 



1 



-G 



--d. 



-GG^'Z' {z) ] + j^Z {z) - G^'g^Al {z) Z {z) 

Z {z) G^ 



(3.9) 



(3.10) 



X 



{X"{x)-g'Al{x)X{x))=0 



where we have used the fact that G^^ = G^^ . In order for this equation to be consistent, we 
must have that: 



X" (x) - g^A' ix) X (x) = -k'X (x) 



where k"^ is a constant. By changing to a diniensionless variable x = ■\/4^jho?x and setting 
X {x) = X (x), equation ( p. 11 ) can be brought to the form: 

X"(x)-^X(x) = -^X(x) , (3.12) 

where k^ = k"^ /2gha'^ . Generically, the solutions to this equation can be written in terms of 
confluent hypergeometric functions, but for the case when k^ is an odd integer, the solutions 
can be written in terms of Hermite functions Hn of order n = (fc^ — 1)/2. We restrict ourselves 
to this case since the Hermite functions decay exponentially for large x, which seems to be 
the natural physical choice. Henceforth in this paper, k^ is taken to be an odd integer, and 
we display some examples of the functions in figure |2|. However, in our present setup, we have 
assumed that the phase of the scalar field remains constant. Since only k^ = 1 corresponds 
to a configuration that preserves its sign, it is the only physical solution for our ansatz. 
As a result, our equation is exactly Schrodinger's equation for a simple harmonic oscillator! 




Figure 2: Some Hermite functions Hn of order n = (fc^ — l)/2, for n = 1, 3, 5 and 101. 



Interestingly, the system confines the x-extent of the condensate in a potential that is exactly 
quadratic. With the constant i?-field passing through the (x, y) plane, we can compute the 
flux of B through a region as: 

$ = BAxAy , (3.13) 

Using the dimensionless variables we introduced earlier, we have: 

, AxAy Ax Ay 



$ = 2/10^- 



(3.14) 



Ax only depends on k'^ since it is entirely determined by the behavior of the Hermite function. 
To estimate Ax, we use that: 

fc2 



(X^) 



c;m2 ~2 



dx \Hn (x)| X = n + 



1 



(3.15) 



and (x) = 0, then our flux becomes: 

! = -#■ (3.18) 

i 2g\l 2 ^ ^ 

where i is the extent in the y-direction, and since it is infinite there, we consider the flux per 
unit y-length. 

Substituting our result of equation ( p.ll ) for the x-dependence back into equation ( p.lOl ), 



the equation of motion for Z (z) becomes: 

^ dJ^/^CZ'{z))+^Z{z)-G''g^A^t{z)Z{z)-2gha'^k^Z{z)G'''' = , (3.17) 



V^^V ' 'J ' L2 

Substituting in the functions, we get: 

It is interesting to note that all a and L dependence in the equation have cancelled. This 
would seem to indicate that the scalar field's behavior does not depend on a and L, but it is 
important to remember that we are working in the perturbative limit where ^ is supposed to 
be small. Earlier, we stated that to ensure the perturbative limit is consistent, we must have 
the quantity aL to be large. Since the only dimensionful parameters are a and L, we take: 

which ensures that the scalar field is always perturbative, and we can now work in terms of 

the dimensionless function Z. 

Near the AdS boundary, equation (p. 18 ) becomes: 



Z" {z) - -Z' {z) + ^Z (z) = , (3.20) 



z ' ' z"^ 



which has the solution: 



limZ(z) = ^iz + *2Z^ , (3.21) 

z—*0 



where ^1 and ^2 are dimensionless constants. Both of these solutions are normalizable, so 
one is not the source of the other. Having chosen m?L'^ = —2 by our choice of potential (see 
equation (|3.lD ), there is not a unique boundary condition at the AdS boundary |2^. ^1 is 
proportional to the vev of the A = 1 operator ((Oi)) and ^2 is proportional to the vev of 
the A = 2 operator ((02)) in the boundary field theory, and only one of these is turned on 
by the boundary condition. This forces us to pick as one of our boundary conditions: 

^1 = or ^2 = . (3.22) 



We choose to only work with the A = 2, since we find the same general qualitative behavior 
for both operators. The exact relationship between the vev and the asymptotic value of Z 
can be calculated using the holographic dictionary (details are shown in Appendix^): 

We can also study the behavior of the solution near the event horizon. We find that there 
are three distinct possible cases. First, if ghk? < 1, then the equation of motion becomes: 

Z" (z) - -^i' (z) + ^Z iz)=0, (3.24) 



where 



9 1 - Ighk"^ 



which has Bessel functions as solutions: 

lim Z (z) = tPiJo (2aV(l-^)) + V'2^0 {2ay/(T^. 

~ ^1+^2 (-7 + -log (a (1-2))) , (3.25) 

\7r vr / 

where 7 is the Euler-Mascheroni constant. Since we want the field to be finite at the event 
horizon, we choose our other boundary condition to be "02 = 0. For ghk^ = 1, the equation 
of motion becomes 

Z"{z)--^Z'iz)+b'^Z{z)=0, (3.26) 



where 



Vg2 



+ . — P^ — >0 ' 



(3_/j2_^2)2 3-/l2-g2 
which has as solutions: 

Ihn Z (z) = iPiJo (6' (1 - z)) + 1P2Y0 [h" (1 - z)) 

-^i + ^2i-l + -\og{h\l-z))\ . (3.27) 

yvr vr J 

Since we want the field to be finite at the event horizon, we would choose the same boundary 
condition as before, i.e. ^^2 = 0. Finally, if ghk^ > 1, we have: 

Z" (z) - r^J' (zq) - ^Z (z) = , (3.28) 

where a^ = 2 g~^f_ ^ > 0, which has as solutions: 



z^l 



hmZ{z) = Vi/o 20^(1-^) +^2i^o [2ay/{l-zq) 



V'i+V'2(-7-log(a(l-^))) . (3.29) 



Again, if we want the field to be finite at the event horizon, we would choose the same 
boundary condition as before, i.e., ip2 = 0. 



4. Numerical Method and Results 

There are several parameters that we can vary in this problem, g, hP', q, h. For the following 
analysis, we take g = 1 for simplicity, but we can expect similar behavior for other values 
of g. In order to solve equation ( p.18 ), we use a shooting method. We impose the following 



initial conditions at the event horizon: 



^w = i. ^'« = |^^ (") 



We find that these initial conditions do not necessarily satisfy equation ( 3.22| ). For a given k'^, 



only for certain values of h and q do we get to satisfy the appropriate boundary conditions. 
Therefore, we fix the value of h, and then scan through possible values of q until the appro- 
priate boundary condition at z = is satisfied, exhibiting our condensate. 

We show in figure y two allowed solutions for the scalar field for a fixed k'^ and h but 
with different q values. Ref. [|l^] argues that only the zero-node solution matters to the phase 
structure, so we consider only these solutions in our subsequent analysis^. It is important to 
emphasize that changing the value of the scalar at the event horizon from the value given in 
equation (J^^ does not change the required value of h and q for the scalar to condense, but 
it does change the value of the vev. Therefore, at particular values of h and q that allow for 
a condensate to form, there is a whole range of allowed vev values for the operator depending 
on its value at the event horizon. For convenience, we define a temperature T (with its 
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Figure 3: Two possible solutions for the scalar field with O2 turned on. The parameters used are 
fc2 = 1 and /i = 0.1. 

corresponding q for which the scalar condenses) to be the temperature at zero magnetic field 
{h = 0) at which we find a solution. This is given by (c/ equation (|2.3D): 



f = J-^. (4.2) 

47r 



*Indeed, we find tliat it is only for tliese zero-node solutions tfiat persist for low enough temperature. 



Using this definition, we can define several dimensionless quantities of interest: 



K4 



y/m 



Att 



L T 




2\/T 



B 

f2 



Att 



2h, -4^ 

L2r2 



47r 



(4.3) 



We present tlie results of our numerical condensate search in terms of these quantities in 
figure 0. 
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Figure 4: The graphs depict the following: the allowed B and T values for O2 to condense, the vev 
for O2, and the allowed p and T values for O2 to condense. 



5. Discussion 

As previously alluded to, since we are working in a limit where the scalar is not back-reacting 
on the Maxwell field, we cannot (as in ref.[^]) track the dependence of the vev on temperature 
all the way to the transition temperature^ Tc- We instead simply exhibit the values of the 
ratios of temperature, charge density, and magnetic field to T that allow a condensate to form 
in our study. This is enough to allow us to study the condensate's spatial behavior, and we 
are firmly below Tc in all that we do.. 



^Note that T is not the transition temperature, but merely a normalization set by the solution at /i = 0. 
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There are several interesting features, as displayed in figure ^. Note that in ref. [|^, the 
scalar condenses at a particular value of T /T^, with T^ oc p^'^ . Similarly, the scalar in our 
setup condenses at a particular value of p/T"^ and B /T"^. There exists a minimum value of 
p/T"^ for the scalar to condense. Beyond these values, condensation occurs but is not believed 
to be stable [|ll[, and the only allowed solution is the trivial one (see the previous section for 
more details). 

At zero magnetic field, the condensate fills the plane and requires the lowest ratio of 
p/T"^. To see that it fills the plane, we recall that the profile in the x-direction for A:^ = 1 is 
given by: 

X(x) = e-T=e-5'^°'^'' . (5.1) 



Therefore, in the limit of /i ^ 0, we have X{x) -^ 1, i.e. there is no x-dependence. As the 
magnetic field is turned on and increases (reading the first and the other two plots from right to 
left as T/T decreases), the value of p/T^ required for the scalar to condense steadily increases 
and the condensate has a finite thickness along the x-direction. At around T/T = 0.25, 
the ratio p/T'^ needed drops rapidly. The condensate also drops rapidly around the same 
interval, which suggests that the magnetic field might be overcoming the forces keeping the 
superconductor together. However, we find that even at T/T = 0, the condensate persists 
with a standard deviation of 1/a/6^q?. 

This result is of particular interest because, in the limit of large magnetic field, i.e. 
a -^ oo, the standard deviation approaches zero. Therefore, the magnetic field, as it grows, 
shrinks the condensate away completely. This is reminiscent of the Meissner effect, where the 
magnetic field expels the condensate. The condensate itself cannot expel the magnetic field, 
as is usually the case, since the scalar cannot back-react on the background magnetic field. 

It is interesting to speculate on what higher k'^ values could mean if for example the 
ansatz for the scalar field is modified to make them physical. As an approximation, we can 
naively proceed with our setup with higher values of k'^. We find that as k'^ increases and the 
x-width of the condensate expands (see figure g), the corresponding B-field associated with 
it is smaller. It is interesting to follow this to large k'^, the "classical" limit of the quantum 
harmonic oscillator controlling the x-profile: The fc^ — > oo limit has the condensate filling the 
plane while the magnetic field B ^ 0. Pleasingly, this is consistent with the limit of small 
and large magnetic field we discussed above^. 
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A. Holographic Dictionary 

The solution of the equation of motion admits two normalizable solutions at the AdS bound- 
ary, which we reproduce here: 

limZ(2) = ^iz + *2^^ , (A.l) 

To proceed with calculating the vev of the A = 2 operator O2 using the holographic dictionary, 
the procedure is to assume that ^1 is the source of the operator [^. To proceed, we first 
write the asymptotic solution of the full scalar field as: 

\mi^!{x,z) = e''^ i(t)Q{x)z + A{x)z'^) , (A.2) 

2:— ►0 

where we have explicitly shown the constant phase of the scalar field. Next, we calculate the 
variation of the on-shell action: 



<J5on~shell = -^ [d^xV^G''d,-^{x,z)6^{x,z) 



2 = 1 



2 = 



= lim — TT / (fx — ^ (z(j)o(x)6(l)o(x) + 2z'^A(x)5(j)o(x) + z'^(j)o(x)5A(x) + 0(z^)) 
2-^0 2^4 J z^ 

At this point, we may worry about the divergence produced by the first term in this expression, 
but we have not included the following counterterm in the expression of our action: 

^counter = -T^V^T / d^^^^i^, Z = 0)^(x, Z = 0) , (A.3) 

4k| L J 

where 7 is not the Euler-Mascheroni constant. Varying this counterterm and including it in 
the on-shell action gives as a final result: 

<55on-shcii = 7r^ d^xL'^a^A{x)6(l)o{x) (A.4) 

Therefore, we have our final result: 



1 OOon-shell L ^2 3 ^/ \ ^ 

-r- ^ , , , = ttL a A(x) = — ■„ 

dp 5(l)o{x) 2k| ^ ' 2kI 



{02{x)) = -;- '^7°7^f^ = f^L2a3^(x) = ^a^m^Xix) , (A.5) 



where d represents the fact that we want to study the operator in terms of unit length in the 
y-direction. If we wish to only consider the overall scale of the operator [i.e., dropping the x 
dependence), we can simply write: 

(O2) = 7^«'*2 . (A.6) 
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